Decrypting the cyclotron effect in graphite using Kerr rotation spectroscopy 
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We measure the far-infrared magneto-optical Kerr rotation and reflectivity spectra in graphite and 
achieve a highly accurate unified microscopic description of all data in a broad range of magnetic 
fields by taking rigorously the c-axis band dispersion and the trigonal warping into account. We find 
that the second- and the forth-order cyclotron harmonics are optically almost as strong as the fun- 
damental resonance even at high fields. They must play, therefore, a major role in magneto-optical 
and magneto-plasmonic applications based on Bernal stacked graphite and multilayer graphene. 
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Owing to its wide spread and technological impor- 
tance, graphite is one of the most studied crystalline 
materials. Magneto-optical spectroscopy was one of the 
key techniques, together with the transport and mag- 
netization measurements, to establish the essentials of 
the unusual band structure of this layered semimetal. A 
small mass, low density and weak scattering of electrons 
and holes give rise to remarkably strong cyclotron res- 
onances [THTT] . Surprisingly, in spite of the large num- 
ber of experiments that appeared after the first report of 
the cyclotron effect [T], the magneto-optical lineshapes 
and intensities are still not satisfactorily explained by 
microscopic calculations, which is in part because of the 
incompleteness of the existing optical data. While the 
measurements where circular polarization served to dis- 
tinguish electrons and holes were typically done with 
only a few radiation wavelengths [2- 5 , most broad- 
band spectroscopic experiments [THSlfTT] provide a highly 
mixed response of electrons and holes. As there ex- 
ist several open questions regarding the optical inten- 
sities in graphene [TMT4] . a quantitative understanding 
of graphite as its three-dimensional counterpart is highly 
desirable. 

In order to combine the advantages of broadband spec- 
troscopy and a sensitivity to the charge carrier type, we 
measure the Kerr rotation angle 6k(&) as a continuous 
function of the photon energy hu, in addition to the more 
conventional reflectivity spectra taken on the same sam- 
ple. The Kerr angle is defined as the rotation of the 
polarization plane upon a normal reflection of linearly 
polarized light from the surface of a sample, when the 
magnetic field B is perpendicular to the surface. The re- 
flectivity is a ratio of the intensities of the reflected and 
the incident beams of unpolarized radiation in the same 
geometry. These quantities are related to the complex 
reflectivity coefficients r±(u>) for the right (+) and left 
(-) circular polarized light: 

_ arg(r_) - arg(r+) _ |r_| 2 + |r+| 2 



and eventually determined, via the usual Fresnel equa- 



tions, by the circular conductivities <j±{lu) — o xx {uj) ± 
ia xy (uj). While the electron and hole cyclotron reso- 
nances are indistinguishable in the reflectivity, they can 
be separated by the Kerr rotation. The latter has the 
same physical origin as the Faraday rotation [Tb\ , but can 
be applied to opaque samples. This set of measurements 
can be regarded as infrared Hall spectroscopy |16j since 
it probes both the longitudinal and the Hall conductivi- 
ties, a xx and a xy , at infrared frequencies as specified in 
the Supplemental Material. 

Figures la and lb show the Kerr angle and the reflec- 
tivity of highly ordered pyrolytic graphite (HOPG) at 
5 K for magnetic fields up to 7 T. The Kerr spectra are 
complicated and highly structured above 1 T; the struc- 
tures displace in energy roughly proportionally to the 
field. They crossover to a simpler shape at lower fields 
before disappearing, as expected, at B — > 0. At higher 
fields the quantized Landau level (LL) transitions can be 
resolved, while at low fields different transitions merge 
and form a broad continuum. The reflectivity demon- 
strates a similar trend, except that it remains finite at 
zero field, and this agrees with previously taken magneto- 
reflectivity spectra at the corresponding photon energies 
and magnetic fields [7J [TT] . 

The extremely small Fermi surface of graphite is elon- 
gated along the K-H line at the edge of the Brillouin zone 
(shown in the inset of Figure 2a), perpendicular to the 
graphene planes. Between K and H, the in-plane bands 
change continuously, the details being strongly depen- 
dent on the relative strength of the different interplane 
hopping parameters (T7J [IS]- However, it is often as- 
sumed in the literature that the observable spectral fea- 
tures originate only from the K-point, where the Fermi 
surface is electron-like, and from the H-point, where it is 
hole-like. An appealing aspect of this simplified view is 
that the band structure for the H-point strongly resem- 
bles the conical band dispersion of monolayer graphene, 
while the bands at the K point disperse parabolically sim- 
ilar to bilayer graphene. The LLs and therefore the cy- 
clotron frequencies at these points are contrasted by their 
dependence on the perpendicular magnetic field: ~ yB 
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FIG. 1. (color online) The Kerr angle (a) and reflectivity (b) 
spectra of highly ordered pyrolytic graphite measured at 5 K 
at various magnetic fields. The curves are offset for clarity; 
the dashed lines indicate at panel (a) and 1 at panel (b). 
Circles are the experimental data, solid lines are the theoret- 
ical fits using the same set of tight-binding parameters at all 
fields. Note that only the Kerr angle was used to adjust the 
band parameters. 

and ~ B respectively [31 [5] • 

We found that it is impossible to reproduce the mea- 
sured optical curves using this two-point approxima- 
tion, i.e. by considering a weighted superposition of 
the theoretical optical spectra of monolayer and bilayer 
graphene |19) . In contrast, taking the entire Fermi sur- 
face into account resulted in excellent fits, shown in Fig- 
ures la and lb as solid lines. Within this approach, the 
total optical conductivity is given by an integral over the 
K-H line [20] 

f H 

cr±(ui) = / a±(k z ,uj)dk z (2) 
Jk 

where a±(k z ,uj) is the partial contribution from all op- 
tically allowed transitions between LLs at a given per- 
pendicular momentum k z . In order to compute the 
LLs and optical matrix elements, we used only the con- 
ventional tight-binding parameters 70 , . . . , 75 and A of 
the well known Slonczewski-Weiss-McClure (SWMcC) 
band model [T7j [TS] and the scattering rate T, which 
we assumed, for the sake of simplicity, to be energy-, 
momentum- and magnetic field-independent. We took 
the SWMcC values from the recent magnetotransport 
measurements |21) and fine-tuned them to fit our Kerr 
spectra at four field values 1,3,5 and 7 T simultaneously, 



thus covering small fields and high fields (almost up to 
the quantum limit 7, 22 ) at the same time. Remarkably, 
the Kerr spectra at other fields as well as all reflectivity 
spectra were reproduced satisfactorily without extra fit- 
ting. The resulting parameters are fairly consistent with 
previous studies [2D [23] ■ Thus in the present case the 
Kerr spectra alone are, in principle, sufficient to deter- 
mine the electronic band parameters. The details of the 
calculations and the fits are given in the Supplemental 
Material. 

One can notice that some of the LL structures, es- 
pecially at low frequencies, are sharper in the experi- 
ment than in the fitting curves. This is likely due to 
our assumption that the scattering rate (T w 2.8 meV 
was found by the least-squares technique) is constant. 
Indeed, reducing T makes the overall fit worse, but it 
improves the match for some spectral features. 

Like it is in bilayer graphene, each LL in graphite (ex- 
cept for the first two) is subdivided into four sublevels, 
corresponding to one of the four 7r-bands, two electron- 
and two hole-like. However, in contrast to graphene, the 
levels in the three-dimensional graphite acquire a k z dis- 
persion [TT] Hi], shown in Figures 2a for 5 T. We 
adopt the notation HH [20] , where the main LL index 
n starts from zero, and the minus sign is used to denote 
the hole sublevels. Additionally, we use the symbol H 
after the index to distinguish the high-energy sublevels. 
The LLs with n = and n — 1 are special as they have 
only one and three sublevels respectively. 

Although the bands and the LLs depend on all the 
SWMcC parameters, the symmetry of the wavefunctions 
and the optical selection rules are determined by the in- 
terplane next-nearest neighbor hopping 73 ss 0.3 eV. It 
is responsible for the so-called trigonal warping [25] [27J , 
which breaks the in-plane radial symmetry of the bands. 
The case of zero warping is much easier to treat analyti- 
cally, since the entire Hamiltonian in a magnetic field is 
factorized into a series of uncoupled 4x4 blocks for each 
LL. If 73 7^ then one has to diagonalize numerically an 
infinite matrix 24J. This case, however, is most inter- 
esting because it comprises nontrivial effects such as the 
magnetic breakdown in graphite |24j and Lifshitz tran- 
sitions in bilayer graphene as a function of strain [28) . 
Furthermore, trigonal warping was shown to have a pro- 
found effect on the charge transport in weakly doped bi- 
layer graphene [29] . 

Without trigonal warping, only the transitions where 
the LL index n changes by ±1 are allowed, due to the 
harmonic-oscillator-like structure of the LL wavefunc- 
tions. Warping mixes LLs separated by 3 times an in- 
teger number (i.e. levels shown by the same color in 
Figure 2a) and makes also the high-order resonances 
An = 3m±l (m = ±1, ±2...) optically active [26]. These 
cyclotron harmonics were indeed observed [fl [21 SI [7J 
and served [21 El [27J to determine 73. However, little is 
known about the absolute optical strength of these tran- 
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FIG. 2. (color online) Panel (a): calculated low-energy Landau-level structure of graphite at 5 T. The inset shows the Brillouin 
zone of graphite, (b)-(k): The theoretical spectra of Re <t+(cj) (right panels) and Re er_(o;) (left panels) at 5 T and 5 K. b, c: 
the total conductivity; (d), (e): the contribution from the fundamental transitions An = ±1; (f), (g): the contribution from 
the high-order transitions An = ±2; (h), (k): the same for An = ±4. The horizontal axis for the left panels is inverted for the 
illustration purpose, since Re a+(—uj) — Re <j_(lj). Conductivities per atomic layer are shown in the units of <to = e 2 /4ft. The 
strongest intraband transitions corresponding to the selection rules n — > n + 1, n — > n + 2 and n — > n + 4 are shown in panel 
(a) by the red, blue and green arrows correspondingly. 



sitions. The very good match between the data and the 
theory allows us to clarify this issue. 

Figures 2b and 2c show the real parts of the total right 
and left-handed optical conductivities as given by the 
model at 5 T and 5 K. The curves are normalized to 
(To = e 2 /4h, since the conductivity of graphite in zero 
field is close to this universal value in a broad energy 
range 30 . The net contributions from the transitions 
satisfying the fundamental selection rule An = ±1 are 
presented separately in Figures 2d and 2e, where also an 
assignment is given to every peak. The most intense res- 
onance at about 13 meV in <r+(w) is due to the electron- 
like LL intraband transitions of the type n — > n + 1. The 
actual number of the transitions involved depends on the 
magnetic field; in particular two of them (1 — > 2 and 
2 — > 3) are activated at 5 T. They stem from the range of 
k z starting from K until approximately the midpoint be- 
tween K and H. This peak is somewhat analogous to the 
cyclotron resonance in usual semiconductors. Addition- 
ally, a series of interband peaks of the type — n — > n + 1 
in <j+(oj) and — n— 1 — > n in cr-{uS) are present [81 |2"UI |2"5] . 
The values of k z involved in these optical transitions span 
the entire Brillouin zone. 

The calculated partial optical conductivities due to 



the high-order resonances An — ±2 and An = ±4 are 
shown in Figures 2f-k. The dominant contributions are 
the electron-like intraband transitions n — > n+2 to <J— (cj) 
and n — > n + 4 to a+ (w) . Although one may expect these 
overtones to be much weaker than the main cyclotron 
resonance, they have, in fact, a comparable intensity and 
form outstanding satellite structures in the total conduc- 
tivity (Figures 2b, c). Qualitatively, one can relate this 
very peculiar observation with the comparable values of 
73 and other interlayer hopping parameters, for exam- 
ple 71 fa 0.4 eV. Interestingly, the optical intensities of 
these two harmonics are about the same, contrary to a 
naive expectation of having a stronger transition for a 
smaller An. This can be attributed to the fact that trig- 
onal warping most strongly mixes LLs separated by 3, 
thus making the harmonics for An = 3 ± 1 the strongest 
and roughly equal in intensity. Our calculations show in- 
deed that the transitions of higher orders (5, 7, etc.) are 
significantly weaker than the first two harmonics. 

Figure 3 presents the evolution of the calculated op- 
tical conductivity spectra as a function of the magnetic 
field. It is evident that the intensities of the transitions 
n — > n + 2 and n — > n + 4 are of the same order of magni- 
tude as the one of the transitions n — > n + 1 up to at least 
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FIG. 3. (color online) The theoretical spectra of Re cr_(w) 
(on the left) and Re <j+(oj) (on the right) at 5 K at the ex- 
perimental values of magnetic field. In the calculations, the 
same parameters were used that give the best fit to the data. 
The curves are shifted for clarity; the zero levels are shown 
by dashed-dotted lines of the same color. The peaks corre- 
sponding to the fundamental cyclotron resonance n — > n + 1 
and the most intense harmonics n — > n + 2 and n — > n + 4 are 
marked by shaded regions. 



7 T and all of them originate from the Drude peak at zero 
field. In fact, the same calculations (not shown) predict 
that the high-order harmonics should remain optically 
visible up to extremely strong fields ~ 100 T. Impor- 
tantly, such a strong effect of higher harmonics must also 
exist in multilayer graphene with a similar (Bernal-type) 
stacking, where trigonal warping is present |19j . 

To summarize, we demonstrated that a classical band 
model, if rigorously applied, describes accurately the cy- 
clotron spectra in graphite in a broad range of magnetic 
fields. This should serve as a solid basis for the optical 
investigation of more subtle phenomena such as the cou- 
pled electron-hole plasma, electron-phonon interactions 
and the spin-orbit coupling, the effects essential for the 
emerging fields of graphite(graphene)-based plasmonics 
and spintronics. We anticipate that measuring the spec- 
troscopic Kerr angle that allows distinguishing the high- 
frequency dynamics of electrons and holes without using 
the circular polarized radiation, will be highly useful in 
the future studies of the topological insulators, interca- 



lated graphites, multilayer graphene and other materi- 
als of fundamental and practical importance. Finally, 
our work suggests that high-order cyclotron resonances 
should also be strong in the optical spectra of multilayer 
graphene up to high fields and therefore should play a 
major role in the related applications. 
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and thank I. Crassee and A. Akrap for critically reading 
the manuscript. 
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EXPERIMENT 

Magneto-optical measurements were done on a large 
piece (about 7x7 mm 2 ) of highly ordered pyrolytic 
graphite (HOPG) of the ZYA grade, the same as used in 
Ref. [Ml The misorientation of the z-axis is smaller than 
0.4° as the X-ray rocking curve showed. The sample was 
mounted in a split-coil superconducting magnet attached 
to a Fourier-transform spectrometer. A mercury infrared 
light source, a Ge-coated mylar beamsplitter and a liq- 
uid helium cooled bolometer detector were used. The 
magnetic field was applied along the z-axis perpendicu- 
lar to the sample and almost parallel to the propagation 
of light. The average angle of incidence is 8°. 

The Kerr rotation angle was determined by placing a 
rotating grid-wire gold polarizer before the sample and 
a fixed one just after the sample (before the detector). 
Without the magnetic field, the signal is at minimum 
when the two polarizers are crossed. In our configuration 
the minimum value was less than 1% of the maximum 
value (when the polarizers are parallel), which allows an 
accurate polarimetric measurement. 

In magnetic field, the minimum shifts by the amount of 
the Kerr rotation. At every field value, we collected the 
reflected spectra for a number of polarizer angles around 
the minimum. The minimum position m j n (w) was then 
determined as a continuous spectrum. The Kerr rotation 
is experimentally determined as 



'K 



(w, B) = (f> mm {u, B) - </> m in(u}, 0), 



(1) 



while the direction of the polarizer rotation was chosen 
to match the physical definition given by Eq. (1) of the 
main text. We verified that 9k(oj,~B) = —8k(uj,B) 
within our experimental accuracy. Therefore, in order to 
reduce the measurement noise, we eventually determined 
the Kerr angle using the relation: 

9 K (u, B) = [4w(w, B) - <t> m in(u, -B)]/2. (2) 

The experimental precision is better than 1 mrad. 

The reflectivity of the sample R{oj, B) was determined 
using a triple-reference method. First, a ratio of the in- 
tensities of a sample (I s ) and a gold mirror (7 r ) was taken 
at a given field B: 



The switching between the sample and the reference was 
done with the aid of a computer controlled translation 
stage. The same ratio p(ui, 0) was taken at zero field. The 
absolute reflectivity was determined using the formula: 



R(u,B) =R(u,0)p(w,B)/p(w,0), 



(4) 



where R(u>, 0) is the zero field reflectivity, measured on 
the same sample in a different cryostat allowing an in-situ 
gold evaporation without displacing the sample [T] . Such 
a procedure eliminates the sensitivity of the detector to 
magnetic field, temporal drifts and other systematic un- 
certainties. 



THEORY 

In deriving the Landau levels in graphite, we essen- 
tially follow the method of Nakao [3] . The main dif- 
ference is that instead of adopting the original SWMcC 
Hamiltonian, we use the tight-binding Hamiltonian with 
equivalent parameters 70, 71, 72, 73, 74, 75 and A, in or- 
der to make a direct connection to the recent literature 
on single- and multilayer graphene The relation be- 
tween the two Hamiltonians and their parameter sets is 
given, for example, in Ref. [5] 

The crystal structure of Bernal stacked graphite is 
shown in Fig. 1, where also the SWMcC hopping pa- 
rameters 7o,...,75 are specified. Due to a strong elec- 
tronic anisotropy, the band dispersion close to the K-H 
line is customarily treated by the tight-binding model in 
the z-direction, and by linear expansion in the in-plane 
momentum (k x ,k y ). For each value of k z one can thus 
define: 

7i,3,4 = 271,3,4 cos(M/2) (5) 
72,5 = 272,5 cos(fc 2 d) (6) 

and write down the 4x4 in-plane Hamiltonian in the 
{Bl, Al, A2, B2} basis: 



H = 



where fc-i 



/ 72 
hvk^ 
hvik+ 



hvk- 
75 + A 

7i 



hv^k 
7i 

75 

hvkj 



A 



hv 3 k + ^ 
hv^k^ 

hvk^ 

72 ) 



(7) 



k 



P( u,B) = I,(u,B)/I r (u,B). 



(3) 



,. v = (3/2)(7 a/fi) « 10 6 m/s is 
the Fermi velocity and where we also introduced £13,4 — 
1)73,4/70. The last SWMcC parameter A stands for 
the on-site energy difference between the non- equivalent 
atomic positions A and B. Here a — 1.43 A is the carbon- 
carbon distance and d = 6.7~A is twice the interlayer 
separation (Fig. 1). 

The velocity operators for the right and the left polar- 
izations can be obtained straightforwardly: 

0i 3 \ 
1 OH _ v 
hdkl ~ V4 



v 

\ V4 v 



(8) 
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order r ,ri,r 2 , 
as a matrix: 



Then the Hamiltonian can be written 
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FIG. 1. (color online) Crystal structure of Bernal graphite 
with the indication of the hopping parameters of the SWMcC 
model. 



1 dH 

hdkZ 



(9) 



In a perpendicular magnetic field B, we make the substi- 
tution 



fc_ ->■ a+^/2eB/h 
k + -> ay / 2eB/h, 



(10) 
(11) 



where a + and a are the ladder operators, acting on har- 
monic oscillator-like wavefunctions "0 n : 



The Hamiltonian in the operator form reads: 



(12) 
(13) 



H 



( 72 ea + €40+ e 3 a \ 
ea 75 + A 71 e 4 a + 
64a 7i 75 + A ea+ 



64a 



(14) 



72 / 



where we introduced the cyclotron energies: 



e = vV2ehB 
£3,4 = e(v3 t4 ,/v). 

The states -0o, •■■ form a complete orthogonal basis 
on one sublattice. As in graphite there are four atoms 
per unit cell, we add the sublattice index m = 1,2,3,4 
and define the following set of basis states: 



I Sim \ 

S 3 m 
\ Sim ) 



ipn- 



(15) 



It is useful to group them in a particular way [2], 
namely: T = {^}, T 1 = {tfj, and T„ = 

{*i,¥jUi,*li_i,¥£_a} for n > 1 and arrange in the 



where: 



tf = 



n=2,3 



Gq 



G n= 2,3... — 




e 3 V2 





/0 e 3V ^T+T\ 





\ / 



The Landau levels are calculated by the diagonaliza- 
tion of the matrix (16). In the case 73 = (no trig- 



onal warping) all G n vanish and the Hamiltonian can 
be factorized into individual blocks H n each of those 
can be solved separately to provide the sublevels. If 
73 7^ then one has to diagonalize the infinite Hamil- 
tonian, which can be further factorized [5] into three in- 
dependent, although still infinite, matrices by separat- 
ing subsets A = {r ,r 3 ,r 6 ,..}, b = {ri,r 4 ,r 7 ,..} and 

C = {T2, T5, Tg, ..}. Numerically, the problem can be 
tackled by truncating an infinite matrix to obtain a fi- 
nite square matrix M x M and making M large enough 
so that increasing it further changes negligibly the en- 
ergies and the eigenf unctions of the low-index Landau 
levels of interest. The bands resulting from subsets A, B 
and C are marked in Fig. 2a of the main text by black, 
red and blue colors respectively. 

The velocity operator in the same representation is 



v+ = v_ = 



( Uo 

Vi u x 

v 2 u 2 

V3 



V 



(17) 
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where 



V x = 



3,4... 



U 




11 



( w 3 \ 



23 ■■ _ 

\0 J 

The circular optical conductivity is calculated using 
the Kubo formula (the spin and the valley degeneracies 
are taken into account): 



8e 3 B 



d**£ \Q\*±\3) 



/(*) ~ f(ej) 



iT 



(18) 



The indices i and j run over all Landau levels at a given 
k z . Here /(e) = (1 + exp (e — ^)/T) _1 is the Fermi dis- 
tribution function. The chemical potential /i depends on 
the magnetic field and can be adjusted using the electri- 
cal neutrality condition [TU] : 



/ " dk z V [fia)- 1/2] =0. 
Jo 



(19) 



The optical longitudinal and Hall conductivities are 
derived from the circular one: 

<7xx(w) = [c+M + a-((jj)]/2 
& X y{u) = [o+M - <r-(u))]/2i. 

In the absence of trigonal warping, only the matrix 
elements between the levels which quantum numbers n^ 
and rij differ by ± 1 are non-zero. If 73 7^ then in 
principle all transitions are allowed except those where 
n changes by 0,±3,±6..., i.e. when both LLs originate 
from the same subset (A, B or C). As our calculation 
shows (see the main text) the transitions with An = ±2 
and ±4 have the strongest optical intensity, while the 
others (for example An = ±5 and ±7) are much weaker. 

In reality one has to limit the summation to a fi- 
nite number of LL transitions. Since our goal is to 
calculate optical properties in a limited spectral range 
< to < Umax, it is natural to keep only the resonances 
below a certain threshold value fl > uj m ax- Obviously, 



the total number of counted transitions then depends on 
the magnetic field. 

The reflectivity and the Kerr angle are given by the 
formulas 



R = 



\r + \ 2 + \r\ 2 



IK = 



arg(r_) - arg(r+) 



(20) 
(21) 



where r± are the complex reflectivities for the circularly 
polarized light. At a normal incidence, they are given by 
the Fresnel equation: 



r± = 



1 + 



where 



(22) 



(23) 



is the in-plane dielectric function. The purpose of the pa- 
rameter is to account for the static dielectric constant 
associated with all higher-energy optical transitions, not 



covered by Eq. ( 18 ). It should not be considered as a free 



parameter, since it is determined by the choice of fl and 
can be fixed using the experimentally obtained dielectric 
function at zero field pQ. 

The integration over k z in practice is done by a sum- 
mation over a finite grid. One has to choose it in such 
a way that making it denser does not change the result 
anymore. The optimal density depends on scattering (r), 
temperature (T) and field (£>). 



FITTING PROCEDURE AND RESULTS 



We fitted the Kerr spectra at 1, 3, 5 and 7 T simulta- 



neously with the Equations (21), (22), (23) and (18), by 



adjusting the SWMcC parameters except 70 that we fixed 
to a value corresponding to v = 1.02 x 10 6 m/s fRef . ITTj) . 
All parameters are assumed to be field-independent. 
The least-squares criterion was used and the Levenberg- 
Marquardt fitting routine was employed [12] , with an an- 
alytical calculation of the derivatives of the output with 
respect to parameters. The RefFIT software [13] was 
used with a dedicated model for graphite. Only a few it- 
erations were needed to achieve a convergence. The best 
values are presented in Table [H] where also the results of 
the fitting of recent magnetotransport measurements [14] 
are shown, as well as the SWMcC parameters obtained 
by fitting the calculated bands using the first-principle 
G W method [TS] . One can see that our parameter values 
agree well with the previous works. The relatively large 
error bars for 74, 75 and A are due to the fact that their 
effects on optical spectra are correlated. By fixing one 
of them to other experiments, or using magneto-optical 
data at higher frequencies and magnetic fields [151 E] 



> 
E 




FIG. 2. (color online) Magnetic field dependence of the chem- 
ical potential. 



TABLE I. SWMcC parameters (in eV) obtained in the present 
work (1st column), in Ref. 1141 (2nd column) and Ref. 1151 (3rd 
column) . 





This work 


Schneider et al. 


1141 Gruneis et al. 15 


7o 


3.16 (fixed) 


3.37 ± 0.02 


3.05 


7i 


0.38 ± 0.01 


0.363 ± 0.05 


0.403 


72 


-0.0089 ± 0.0003 


-0.0121 ± 0.0005 


-0.0125 


73 


0.297 ± 0.005 


0.31 ± 0.05 


0.274 


14 


-0.15 ± 0.1 


-0.07 ± 0.01 


-0.143 


75 


0.005 ± 0.03 


0.025 ± 0.005 


0.015 


A 


0.027 ± 0.05 


0.024 ± 0.002 


0.05 


r 


0.0028 ± 0.0002 







TABLE II. The used computational settings for every field 
value as explained in the text. 



B (T) 0.5 


1 


1.5 


2 


3 4 5 6 
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M 180 

N z 50 
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FIG. 3. (color online) Real parts of a xx (u>) (a) and a xy {uj) 
(b) that correspond to a±(cj) of Fig. 3 of the main text. 



we can likely reduce the error bars. However, the main 
goal of the present work is to demonstrate that the Kerr 
spectra are fully consistent with the SWMcC model. 

Fig. [2] shows the dependence of the chemical potential 
on magnetic field. The non-monotonic dependence is due 
to the crossing of the Fermi level by the LLs. Almost 
identical curve is obtained in Ref. [21 

In the calculations, both f2 and €00 were the same at 
all fields (0.2 eV and 15 respectively). The size of the 
truncated matrix M and the number of k z grid points 
N z was optimized for every field in order to speed up 
the calculations while keeping a sufficient accuracy. The 
actual values used are given in Table II. 

For completeness, in Fig. [3] we present the longitudi- 
nal and the Hall conductivities obtained from the same 
calculation (Fig. 3 of the main text). 
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